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Abstract 

We report the experimental verification of noise-enhanced logic behaviour in an electronic analog 
of a synthetic genetic network, composed of two repressors and two constructive promoters. We 
observe good agreement between circuit measurements and numerical prediction, with the circuit 
allowing for robust logic operations in an optimal window of noise. Namely, the input-output char- 
acteristics of a logic gate is reproduced faithfully under moderate noise, which is a manifestation of 
the phenomenon known as Logical Stochastic Resonance. Interestingly, the two dynamical variables 
in the system yield complementary logic behaviour simultaneously, indicating strong potential for 
parallel processing. 
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I. INTRODUCTION 



Realization of logic functions in different physical systems is one of the key questions 
that commands widespread research interest in sciences and engineering. Universal general- 
purpose computing devices can be constructed entirely from NOR/NAND logic gates It 
is particularly interesting to investigate if systems of biological relevance can also yield logic 
outputs consistent with the truth tables of different logic functions (see Table [T]). 

A new idea in this direction uses the interplay between noise and nonlinearity construc- 
tively to enhance the robustness of logic operations. Namely, in an optimal window of noise, 
the input-output characteristics of a logic gate is reproduced faithfully. This phenomenon 
is termed Logical Stochastic Resonance |2|, |3|. 

In this work we investigate the possibility of obtaining such reliable logic outputs exper- 
imentally, in an electronic analog of a noisy biological system of considerable interest: a 
synthetic genetic network {4]. We demonstrate the pivotal role of noise in the optimization 
of the logic performance in this circuit. Further, we show the capacity for parallel processing 
... w.. co.p,e.e.a. ,o,c ou.pu. oMai„ed .^.^o^ ^ . 

Since understanding the intercelluar processes in a network of interacting biomolecules 
is difficult, an alternative approach has been started recently jol, to design artificial genetic 
networks to derive desired functional behaviors. One important early design is a clock using 
three genes inhibiting each other in a cyclic order. Each gene produces a repressor to the 
subsequent gene. The genes each produce their own mRNA, which translate the repressor 
protein. Taking into account the standard chemical kinetics for production, degradation 
and inhibition, a dynamical system model was proposed [3] where the repressor-protein 
concentrations and mRNA concentrations were expressed as dynamical variables. 

Here we use two such repressors and constructive promoters as our model system for 
implementing logic functions. First we describe the synthetic gene network model below, 
and define what constitutes logic inputs and logic outputs in this system. We then go on to 
present the electronic analog of the system followed by a comparison of numerical simulation 
and experimental measurement. 
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TABLE I: Relationship between the two inputs and the output of the fundamental OR, AND, NOR 
and NAND logic operations. Note that the four distinct possible input sets (0,0), (0, 1), (1,0) and 
(U, ..... c„„d.t„™ a. (0. 1) and (, 0, Note t>„. any l„..ca> c...t 

can be constructed by combining the NOR (or the NAND) gates [1|. 

II. SYNTHETIC GENETIC NETWORK MODEL AND LOGIC OPERATION 



We consider the model of two genes inhibiting each other. The concentrations of the two 
expressed proteins are x and y, and their rates of change are: 
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where [32 are the rate of decay of each expressed protein and n is the Hill coefficient. 
The «!, a2 describe the maximum transcription rate in absence of inhibition and they are 
used here as tunable parameters. In the original model gi and g2 represent leakiness of the 
promoter. However we use them as constant bias. The additive noise has strength D and 
r] and ^ are chosen from unit variance zero mean Gaussian distributions. Ji and I2 are two 
low amplitude inputs. 

The system above has two stable configurations in the x ~ y plane: one state has a high 
value of X {xu) and a low value of y (yi); the other state has a low value of x (xi) and a 
high value of y {yu)- That is, the two dimensional potential underlying this system has two 
wells, {xu,yi) and {xi,yu), in the x — y space. Varying the parameters changes the depth 
and position of these wells. 

Encoding Inputs: Here the low amplitude input signal is / = Ji + I2, with /1//2 equal 
to Ion {Ion > 0) if the logic input is 1, and /1//2 being if the logic input is 0. So we have: 

(i) Ii + I2 = corresponds to logic input set (0, 0) 

(ii) Ii + I2 = Ion corresponds to logic input sets (0, 1)/(1, 0) 



(iii) Ii + I2 = 2I0N corresponds to logic input set (1, 1) 



Output : 

The outputs of the system are determined by the level of the dynamical variables x{t) 
and y{t). For instance the output can be considered a logical 1 if the state is at the high 
level, and logical if it is at the lower level. That is: 

(i) If X < X*, then Logic Output is 

(ii) If X > X*, then Logic Output is 1 

Here x* is the output determination threshold that lies between the two states, e.g., at 
the position of the barrier between the wells. The results presented here are not sensitive to 
the specific value of x*. 

Specifically, in this work, we consider the logic output to be 1 when the state is close to 
the upper well, and when the state is close to the lower well, for both x and y variables. 
So when the system switches wells, the output is "flipped" or "toggled" . 



III. CIRCUIT REALIZATION 

t — I 

The circuit of a single inhibitory gene [8l] is shown in Fig. [H The transistor current 
represents the rate of gene expression and the voltage Vout represents the concentration of 
expressed protein. The Vdh adjusts the binding constant of the repressor to the gene's DNA. 
The Hill function inhibition in Eqs. 1-2 is accounted for by the dependence of the transistor 
current on input voltage Vm which represents the repressor concentration. The synthetic 
genetic network shown in Fig. [2] is comprised of two individual gene circuits connected in a 
loop, each inhibiting the other. For the model in Eqs. 1-2, the encoding inputs Ji and I2 
add to production of y which is accounted for in Fig. |2] by the two logic driven transistors 
sourcing current to Vy. 

The circuit equations are obtained by applying Kirchoff 's laws to and Vy, the voltages 

Multiplying both equations by Ry results in 



across the capacitors in Fig. [2] as in Ref. 
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equations for Vx and Vy] 



RyC 




+ Ryit 
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(4) 



where it is the transistor current in Fig. [T] for each gene, and ii and 22 are the logic train 
transistor currents in Fig. |2l The variables (x, ?/) in Eqs. 1-2 are related to voltages Vx and 



where Vth corresponds to the repressor's equilibrium binding constant. The connections 
between model parameters {ai,(3i,n) and circuit parameters are found by relating circuit 
Eqs. 3-4 to the model Eqs. 1-2 and by adjusting the dependence of the transistor current it 
on Vin to match the Hill function inhibition (see Appendix for details). 

A noise generation circuit shown in Fig. [3] based on breakdown of a reverse biased base- 
emitter junction produces noise with zero mean and variable amplitude. We use a well 
regulated supply for the noise circuit to avoid adding AC frequencies into the noise. Two 
of these noise circuits are used to supply noisy voltages to each individual gene at locations 
indicated in Fig. [21 The amplitude D of the gaussian noise in the model is related to the 
rms value of the noise voltage hj D = Vrms/Vth- 

IV. RESULTS AND DISCUSSION 

FigureHlshows measurements from the circuit using n = 2.4, a = 1.78, f3i = 0.90, (32 = 1, 
iHigh = 0.062 for different strength of noise. It is apparent that for noise levels within an 
optimal range (Fig. 4c) the circuit indeed performs the logic AND/NAND function. These 
measurements agree well with the numerical simulations of Eqs. 1-2 shown in Fig. [5] using 
the Euler-Maruyama method with a step size of At = 0.025. 

One can also reconfigure the system to another set of logic functions, namely the 
fundamental OR/NOR logic, by simply changing the value of g2- For instance when 
g2 = 0.06 — 0.07 (with all other parameters unchanged), one gets a clear OR and the 
complementary NOR response from the system. Inclusion of g2 in the circuit is achieved by 
including a current sourcing transistor in the same way as the encoding signal transistors. 



Vyhy. 
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but with the logic train input being grounded. The same analysis as for the encoding values 
and currents (see Appendix) finds that g2 = 0.062. Figures |6] and [7] demonstrate this capacity 
of reconfiguration in simulations and circuit experiments respectively. 

In summary, our results show that the dynamics of the two variables x and y with 
g2 = 0, mirror AND and the complementary NAND gate characteristics. Further, when 
92 7^ 0, we obtain a clearly defined OR/NOR gate. Since x is low when y is high and 
vice-versa, the dynamics of the two variables always yield complementary logical outputs, 
simultaneously. That is, if x{t) operates as NAND/NOR, y(t) will give AND/OR. This is 
significant, as it allows the system to yield two logic operations in parallel, thereby doubling 
the computational speed. 

So our results here extend the scope and indicate the generality of the recently observed 
phenomena of Logical Stochstic Resonance (LSR) through experimental verifications. Fur- 
ther, these observations may provide an understanding of the information processing ca- 
pacity of synthetic genetic networks, with noise aiding logic patterns. It also may have 
potential applications in the design of biologically inspired gates with added capacity of 
reconfigurability of logic operations and parallel processing. 

Thus the results presented in this work suggest new directions in biomolecular comput- 
ing, and indicate how robust computation may be occuring at the scale of regulatory and 
signalling pathways in individual cells and even within few biomolecules. Design and engi- 
neering of such biologically inspired computing systems not only present new paradigms of 
computation, but can also potentially enhance our ability to study and control biological 
systems 
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Appendix 

We give a brief description of the analysis of the single gene circuit in Fig. [H with 
more details in Ref. ^. The op-amp Ul is configured as a subtraction amplifier with gain 
Gi = = 1.1. A replication of the Hill function behavior is achieved by allowing saturation 
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of the output of the op-amp U2 and by having different unsaturated gains (7+2 and G'_2 for 
Vin > Vcth and Vm < Vdh, respectively, due to the diodes in the feedback for U2. G-2 is 
the gain of U2 when its output goes negative, in which case the diodes are not conducting, 
and therefore (j_2 = ff = 3.3. G^2 is a diminishing gain when the output of U2 becomes 
increasingly positive causing the diodes to go into conduction. An increasing repressor 
concentration corresponds to Vm surpassing Vdh which causes the unsaturated output at U2 
to change from a negative voltage of GiG-2(Vi„ — Vct^) to a positive voltage G'iG'+2(^n — Kth)- 
The increasing voltage at the output of U2 turns the transistor off {it — ?■ 0) which corresponds 
to inhibition of protein expression. Transistor current it is maximum when the output of U2 
is at its lower saturation value of V^sat = —3.5 V (for the LF412 op-amp supplied by ±5 V) 
causing about a 0.65 V drop across the 220Q and therefore the maximum current imax = 3 
mA. 

A comparison of Eqs. 1-2 with Eqs. 3-4 yields js]: «i = 03 = A = = li = 0-90' 

/32 = 1. The voltage imaxR is (3 mA)(0.45 kQ) = 1.35 V, where 0.45 kQ comes from the 
Ry = 470n being nearly in parallel with the resistance (10 kQ) at the input to the subtraction 
amplifier Ul for gene-x. The characteristic time here is RyC = (470) (10^^) = 4.7 /xs, and 
r = is the dimensionless time. 

The relation for Hill coefficient n is found by adjusting the dependence of it on Vin to 
match the slope of the Hill function 1/(1 -|- x") at a; = 1 resulting in the relation 8|: 

fVthGiG-2 _ ^ /j-x 
/(5 - V^,at) - 0.6 ~ 4 ^ ^ 

In Fig. [1] the voltage divider fraction / = 0.4/2.6 = 0.154 and V-sat = —3.5 V. Using 

« = imaxR/Vth in Eq. 5 yields the parameter relation na = 1.17GiG-2- To find the relation 

between Vcth and Vth we note that the Hill function equals 0.5 when x = 1. Therefore it 

must be half its maximum value when Vm = Vth which gives 

It f{5-GiG.2{Vth- Vcth)) -0.55 



imax f{5 - V-sat) - 0.6 



0.5. 



Solving gives V^h ~ Vth + q_^q Here we use n = 2.4 and a = 1.78 which givesi G\G— 
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(2.4 X 1.78)/1.17 = 3.65, satisfied by Gi = 1.1 and G_2 = 3.3; Vh = 1.35/1.78 = 0.76 V; 
and Vcth = 0.76 + 1/3.65 = 1.03 V. 

The encoding signals /i, I2 in Eqs. 1-2 are related to transistor currents by: 

_ Ryil.2 

^1,2 - —TT — • 

Vth 



The encoding currents ii^2 in the transistors in Fig. |2]take on two possible values depending 
on whether their logic train input is high or low. When the input is high (> 4V) the transistor 
is off so the current is zero. When the input is zero, there is IV across the Rj = 10 kQ 
connected to the emitter of each pnp transistor creating current ii 2 = (1 V)/Ri. Therefore 
the high level for each encoding input is {Ii^2)High = {Ry x ^y)/{Riyth)) e.g. for a = 1.78 
then Vth = 0.76 V and the high level is 470/(10kx0.76)= 0.062. 
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FIG. 1: Circuit for single gene. Inhibitory input at Vin- Expressed protein concentration is 
represented by Vgut- Ri = 470 ^ for gene-y, 520 for gene-x. Dual op-amp is LF412 supplied by 
-I-/-5 V. The pnp transistor is 2N3906. The input noise has a mean of V (gnd) and controllable 
standard deviation. 
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FIG. 2: Circuit for synthetic genetic network. Encoding inputs are to 5 V pulse trains. The x 
and y gene circuits are shown in Fig. [H Each noise input has its own noise circuit shown in Fig. [3l 
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FIG. 3: Noise circuit and its connection to resistor of gene circuit. LF412 dual op-amps supplied 
from +/-12 V regulators. 
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FIG. 4: Experimental measurements from circuit (cf. Eqs. 3-4) with different noise amplitudes. 
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(a)D = 0.02. 



{h)D = 0.05. 
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FIG. 5: Simulations of Eqns. 1-2, with different noise amplitudes. 
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FIG. 6: Simulation of Eqns. 1-2 with g2 = 0.07 showing OR/NOR logic response. Here noise level 
isD = 0.11. 
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FIG. 7: Experimental measurements from circuit with §2 = 0.062 displaying OR/NOR character- 
istics (cf. Eqs. 3-4). Here noise level is D = 0.10 (i.e. 75 mV rms). 
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